Particle density in zero temperature symmetry restoring phase 
transitions in four-fermion interaction models* 
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Abstract 

By means of critical behaviors of the dynamical fermion mass in four-fermion interaction models, 
we have shown by explicit calculations that when T — the particle density will have a discontinuous 
jumping across the critical chemical potential fj, c in 2D and 3D Gross-Neveu (GN) model and these 
physically explain the first order feature of corresponding symmetry restoring phase transitions. For 
second order phase transitions in 3D GN model when T — » and in 4D Nambu-Jona-Lasinio (NJL) 
model when T = 0, it has been proven that the particle density itself will be continuous across /x c 
but its derivative over the chemical potential /i will have a discontinuous jumping. The results give 
a physical explanation of implications of the tricritical point (T,fj,) = (0, fJ. c ) in 3D GN model. The 
discussions also show effectiveness of the critical analysis approach of phase transitions. 

PACS numbers: ll.10.Wx; 12.40.-y;11.30.Qc 
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1 Introduction 

Four-fermion interaction models E] are good laboratories to research dynamical symmetry restoring 
phase transitions at high temperature and high density It has been proven that when temperature 

T goes to zero, the feature of the phase transitions in these models strongly depend on dimension D of 
time-space jB] It is shown that that the symmetry restoring phase transitions at high fermion 
chemical potential /i and T — are first order when D = 2 and D = 3 however second order when D = 4 
and the momentum cutoff of the loop integrations is large enough. The above conclusions can come 
from analyses of critical behaviors of the dynamical fermion mass m as the order parameter of symmetry 
breaking. The critical behaviors of the order parameter could show the essential characteristics, however, 
they do not directly explain the physical realization of the above phase transitions. The latter is just 
what we want to explore in this paper. Since particle number density is a thermodynamical quantity 
directly relative to the chemical potential /i, we will calculate it in variant cases and combine the results 
with the critical behaviors of the order parameter obtained in Ref. El so as to further expound physical 
characteristics of the above first and second order phase transitions. 

The total particle number density no(T, fJ.-,m) in D dimensional time-space should be a difference 
between fermion and antifermion density and at finite T and finite /i it can be expressed by [S] 
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,y = /3to, a = /x/to, (2) 



where to is the fermion mass and in a four-fermion interaction model with spontaneous symmetry break- 
ing, it will be replaced by the dynamical fermion mass to = m(T, /x) caused by the bilinear fermion 
condensates 0. We will discuss the cases of D = 2, 3 and 4 respectively in Sects. EHH]andE] then come 
to our conclusions in Sect. [5] 



2 D=2 

First let us consider the case of to = 0. We may find out from Eq.(l) that 



n 2 (T, /x, to = 0) 
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Eq. (3) represents the density of massless fermions or the particle density in the case without spontaneous 
Next turn t 

n 2 (T, /x, to 



symmetry breaking. Next turn to the case of to 7^ 0. Let z = y/x 2 /y 2 + 1, we may get from Eq.(2) that 

fi ^ 'Vz 2 ^! _ey( z ~ a ) + 1 
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We will concern ourselves only with the results in the limit T — > (or y — > 00). It is easy to see that 
when y — > 00, the second term in Eq.(4) becomes zero and the non-zero contribution comes only from 
the first term when a > 1, owing to the interval of the integral variable z is from 1 to 00. Thus we may 
obtain from Eq.(4) that 



n 2 (T = 0,/x, to) 




when a < 1 
when a > 1 



(5) 



Eq.(5) is applicable to the general case of free fermions with to 7^ 0. Obviously, ri2(T = 0, /i, to) usually 
varies continuously when fi crosses over m. However, in a four-fermion interaction model with spontaneous 
symmetry breaking and restoration, there could be different case. It has been proven that 0, in a 2D 
GN model, when T = the order parameter m will varies in the following form: 



m(T = 0, n) = to(O)0[to(O) - fi\ 



(6) 



where to(0) is the dynamical fermion mass at T = fi = and just the critical chemical potential [i c at 
T = 0. Substituting Eq.(6) into Eq.(5), we will reach the conclusion that 



n 2 (T = 0, /x, to) 



when fi < to(0) 
/x/7r when /x > m(0) 



(7) 



Eq.(7) indicates that the jumping of the order parameter to from m(0) to across over /x c = m(0) leads to 
jumping of the thermodynamical particle density n 2 (T = 0, /x, to) from to /x/to, i.e. from the zero value 
with spontaneous symmetry breaking to the density of free fermion with to = after the symmetries 
are restored. This clearly shows physical characteristic of the first order phase transition at T = and 
/x = fj, c in 2D GN model. 

3 D=3 



In this case, we will first obtain from Eq.(2) that 
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Eq.(9) will lead to the T ^ limit of n 3 (T, /i, to) 

n 3 (T = 0,/x,m) = ^-9(fi-m)(fi 2 -m 2 ). (10) 

47T 

Assuming that y ~ m/T ^ 1 (e.g. near a critical point of a second order phase transition), then we can 
expand n 3 (T, /x, m) expressed by Eq.(8) in power of y, keeping only the terms up to the order of y 2 , and 
obtain 

Tfl 2 G r — 1 

n 3 (T,n,m) ~ n 3 (T,/j,,m = 0) - — — - (11) 
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where 



n 3 (T,/x, m = 0) = — — — sinh(fc^) 
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Noting that, different from the case of D = 2, the total particle density is temperature-dependent and 
when T = we have 

n 3 (T = 0,/x,m = 0) =n 2 /4n. (13) 

Eq.(10) is a general expression for free fermions with to ^ at T = 0. However, for 3D GN model with 
symmetry restoring phase transition when T — 0, we have the order parameter to whose critical behavior 
is similarly expressed by Eq.(6) fT]- Hence we obtain 

_ x f when u < m(0) ... 
v ^ ' [_ xt 2 /47r when /x > to(0) v 7 

i.e. the zero temperature particle density n 3 (T = 0, /x, m) will jump from with spontaneous symmetry 
breaking to /x 2 /4-7T of massless fermions (symmetries being restored) when /x crosses over the critical point 
to(0) of the phase transition. Such jumping of the particle density definitely indicates the first order 
feature of the phase transition. 

On the other hand, when T ^ 0, near a critical chemical potential /i c , the order parameter to will 
continuously vary in the form [Jj 

to 2 ~ 2Tsinh ^ (/x c - /x), when T^0. (15) 

Substituting Eq.(15) into Eq.(ll) we will obtain 

w— ;> ~ \ sinh fp f/x c — /x) when/x~/x c 

^7T e + 1 i v ^ r~ r~ (16) 
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and 



dn 3 (T,n,m) ^ dn 3 (T, fj,,m = 0) f ^ ^ - j- sinh ^ when ix ~ /x c 
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sases 

from the value less than n 3 (T, /x,to = 0) up to n 3 (T, fi,m = 0), however the derivative dn 3 (T, fj,,m)/d[i 



dpL dfi { when /x > /x c ' 

It is seen from Eqs.(16) and (17) that at T 7^ when /x crosses over /x c , n 3 (T, /x, to) continuously increases 
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has a discontinuous jumping at /i = /i c . This shows the second order of the phase transition. We empha- 
size that this feature can be maintained when T is finite and very small as long as fi c — fJ- is also very 
small so that m 2 /T 2 <C 1 is kept. The limiting point will be (T = 0, /x = m(Q)) which is the same as 
the first order phase transition point when T = 0. The above discussions give a physical explanation of 
implications of the tricritical point (T, //) = (0, m(0)) in 3-D GN model: The particle density will jump 
from to /j 2 /47r when the point (0, 777(0)) is approached along the (T = 0, /j) axis but will continuously 
increase from a lower value to /j, 2 /Air when the point is approached along the T — fx critical curve of the 
second order phase transitions. 



4 D=4 

For AD NJL model, when m = we will have 

1 f°° 

n 4 (T, n, m = 0) = ^2 J d PP 2 
The equivalent differential equation is 
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with the boundary conditions 



n 4 (T,n = 0,0) = 0, 



dn 4 (T,fi,0) 
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The solution of Eq.(19) satisfying the conditions in Eq.(20) is that 

n 4 (T,^0) = — + -n = —. 
When m ^ and T — > (y — * oo), we may obtain from Eq.(2) that 

3 /* oo 
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Eq.(22) is a general result of free fermions and usually does not involve any spontaneous symmetry 
breaking and restoration. Now consider the case in which symmetry restoring phase transition is assumed 
to occur. Thus we should identify m in Eq.(22) with the dynamical fermion mass. It has been proven 
that in a AD NJL model, the dynamical fermion mass to the following behavior [S] 



to = m(0), 
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where A is the 4 dimensional Euclidean momentum cutoff and the critical chemical potential /j c o is defined 
by 



Mc0 = ^ 2 (0)ln 



A 2 



We will assume that 
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so that all the phase transitions will be second order [S]. Combining Eq.(22) with Eq.(23) we may find 
that when fi < m(0), n^T — 0,fi,m) — and when m(0) < [i < /i c o, ni{T = 0,/^,m) will increase in 
the form /i 3 (l — m 2 //z 2 ) 3/ ' 2 /67r 2 (noting that dm/dfj, < in this region). As soon as \i is arriving at /i c 
where symmetries will be restored, m will reduce to zero and (T = 0,/i, m) becomes /i 3 /67r 2 which is 
just the number density of massless fermions given by Eq.(21). It should be emphasized that, since m 
reduces to zero at [i c continuously, n^T = 0, /i, m) changes from into /i 3 /67r 2 also in continuous form 
when /x varies from m(0) to fi c and this is just the characteristic of a second order phase transition. The 
above result can be outlined as 



n^{T — 0, /i, m) 




when m(0) < /i < fi c o 
when n > /i c q 



(26) 



We can also prove that the derivative of n^(T = 0, /i, m) over \i will be discontinuous when /i crosses over 
/i c o. To this end, we will use the critical behavior of the squared order parameter m : 



M 2 ) 
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(27) 



where 6(/i c0 ) = 1/ ^ln - 2) > ^ ow ing to Eq.(25). Substituting Eq.(27) into Eq.(26), we get 



n 4 {T = 0,[i,m)= ' 67r 
and furthermore, 
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Eq.(29) indicates that dui{T = 0, fj,,m/d^, has a jumping over /1 = [i c o and is discontinuous. The above 
results verifies that the discussed symmetry restoring phase transition in the AD NJL model is second 
order indeed. 



5 Conclusions 

By combining the calculated fermion number densities as a thermodynamical quantity with the known 
critical behaviors of the dynamical fermion mass in the four-fermion interaction models, we have shown 
that when temperature T = 0, across over the critical chemical potential /i c , the particle number density 
will have a discontinuous jumping in 2D and 3D GN model and these physically explain the first order 
feature of corresponding zero temperature symmetry restoring phase transitions. The second order phase 
transitions in 3D GN model when T ^ and in AD NJL model when T = and /i 2 < A 2 /4e are 
illustrated by the facts that the densities themselves are continuous but their derivatives over chemical 
potential [i have a discontinuous jumping. These results clearly show physical difference between the first 
order and second order symmetry restoring phase transition in these models, especially give a physical 
explanation of implications of the tricritical point in the 3D GN model. The whole discussions also 
indicate the full effectiveness of critical analysis of the order parameter for phase transition problem. 
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